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BRIEF 


The  theoretical  response  for  the  application  of  square  wave  voltammetry 
at  a dropping  mercury  electrode  is  presented,  with  particular  reference  to 
points  of  current  measurement. 


ABSTRACT 

The  theoretical  aspects  of  square  wave  voltammetry  at  the  droppinq 
mercury  electrode  are  presented.  The  technique  involves  scanning  the 
entire  potential  range  of  interest  on  a single  drop  of  a dme.  Asymmetries 
in  the  waveform  as  well  as  variations  in  current  measurement  parameters  are 
discussed.  Indications  are  that  previous  uses  of  the  waveform  may  not  have 
utilized  all  its  capabilities. 
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For  the  past  several  years,  we  have  been  concerned  with  development 
and  application  of  pulse  polarography  and  voltamnetry  (1).  The  advantages 
of  these  pulse  techniques  reside  chiefly  in  their  discrimination  against 
double  layer  charging  currents,  but,  as  they  are  usually  applied,  they  are 
relatively  slow:  only  one  measurement  is  made  per  drop  of  the  dme  and  a 

ratio  of  delay  time  to  pulse  tim?  of  10:1  or  greater  is  normally  employed. 

Recently  a modification  to  the  PAR  174  was  introduced  (2)  which  allows 
it  to  do  a fast  sweep  mode  of  differential  pulse  polarography.  One  of  the 
problems  inherent  with  this  modification  arises  from  the  linear  ramp  used 
as  the  base  potential  sweeping  function.  The  required  high  slope  of  the 
linear  ramp  reintroduces  charging  current  background  that  the  differential 
pulse  technique  was  intended  to  discriminate  against.  This  fast  sweep 
technique  also  suffers  from  resolution  problems  at  sweep  rates  greater  than 
200  mV/sec  in  that  insufficient  points  are  *”)ken  to  define  the  peaks. 

We  now  turn  to  the  technian*:  ot  squa»e  wave  voltammetry  at  the  dme 
which  retains  a high  discrimination  against  double  layer  charging  current 
but  which  allows  measurement  of  a complete  potential  range  during  the  life 
d?  a single  drop  at  the  dme.  This  technique  was  reported  by  Barker  (3)  as 
long  ago  as  1957,  but  has  sine,'.*  seen  little  application.  Ramaley  and  Krause 
have  developed  the  theory  (4)  and  have  presented  some  application  (5)  of 
square  wave  voltammetry  .t  the  hanging  mercury  drop  electrode. 

The  discussions  cf  Ramaley  and  Krause  were  limited  to  small  step  heights 
(and  consequent  slow  scan  rates).  In  order  to  carry  out  experiments  within 
the  life  of  a single  drop  of  a DME,  faster  scan  rates  (>100  mV/sec)  are 
required.  Faster  scan  rates  imply  relatively  larger  step  heights  for  the 
underlying  staircase  wave-form.  While  the  equations  presented  below  are 
/alid  for  all  values  of  experimental  parameters,  our  calculations  concentrate 
on  those  values  which  are  useful  in  our  particular  case. 
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Our  fundamental  potential  wave-form  differs  superficially  from  that  of 


Ramaley  and  Krause  in  that  the  phase  of  the  square  wave  is  shifted  by  it 


relative  to  the  staircase.  Our  treatment  also  includes  the  possibility  of 


some  asymmetries  in  wave-form  and  measurement  times  not  considered  by 


Ramaley  and  Krause. 


We  develop  our  treatment  for  a stationary  planar  electrode;  it  is 


shown  in  a later  paper  (6)  that  the  effects  of  drop  growth  can  be  well 


compensated  by  an  area  normalization  of  the  observed  currents. 


The  general  square  wave-form,  shown  in  Figure  1,  can  be  thought  of  as 


a pair  of  pulses  superimposed  on  a staircase  wave-form.  One  of  the  pulses 


is  coincident  with  the  step  advance  of  the  staircase,  while  the  other  is 


in  the  opposite  sense  and  lasts  until  the  end  of  the  potential  step.  In 


the  work  of  Ramaley  (4,5),  the  first  pulse  is  opposite  in  sense  to  the 


staircase  direction  while  in  our  case  the  pulse  coincident  with  staircase 


advance  is  in  the  ascending  staircase  scan  direction.  In  the  work  presented 


here,  the  two  pulses  are  of  equal  amplitude  (with  respect  to  the  basic 


staircase)  but  need  not  be  of  equal  duration. 


The  nomenclature  we  use  is  adapted  from  that  of  staircase  voltammetry 


(7)  and  differs  somewhat  from  that  of  Barker  (8)  and  Ramaley  (4,5).  The 


basic  time  unit,  t,  corresponds  to  the  step  width  in  staircase  voltammetry 


and  AE  is  the  step  height.  For  the  square  wave-form,  some  additional  para- 


meters must  be  introduced  as  shown  in  Figure  1.  The  square  wave  amplitude 


is  E$w  and  corresponds  to  the  magnitude  of  each  of  the  two  potential  pulses 


(with  respect  to  the  staircase  potential).  The  fraction  of  each  period 


during  which  the  potential  is  "up"  is  a;  the  length  of  the  forward  pulse 


therefore  is  ot.  The  current  is  measured  twice  during  each  staircase  period; 


at  time  t during  the  forward  pulse  and  at  time  p2  t during  the  reverse 
pulse.  The  coefficient  p^,  p2,  and  a are  all  fractions  less  than  unity  and 
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p-j  < o < pg.  The  treatment  of  Ramaley  and  Krause  was  restricted  to  the 
symnetrlc  case  in  which  a = 0.5  and  p^  = p2 "0i  The  usual  output  is  the 
difference  in  current  measurement  at  p^  t and  p2  t during  the  same  stair- 
case period,  which  we  call  the  forward  difference  and  designate  Al  , but 
other  differencing  schemes  are  possible. 

Assuming  semi-infinite  linear  diffusion  and  a reversible  electrode 
reaction  with  only  the  oxidized  form  initially  present  in  solution,  the 
equation  for  the  currents  can  be  obtained  by  straight-forward  application 
of  the  superposition  theorem  (9)  or  by  generalization  and  modification  of 
the  arbitrary  wave-form  equation  presented  by  Rifkin  and  Evans  (10). 

During  the  ith  period  of  the  staircase,  the  current  measured  at  p^  t 
during  the  forward  pulse  is 


nFAD^  ( C2( j-1  )-C-|  ( j)  Cgtm-l  )-C]  (m)  C^mJ-C^m) 

i 1 = — ■ \ : + 1 ; + r 


Uj) 


r ITT 


I 


p *s 

1 


I 


m=o 

and  the  current  at  p2t  during  the  reverse  pulse  is 


( j-m+p^o)' 


Ir(j)  - 


nFAD, 

/iTT 


J-  j (^(m-D-C^m)  C1(m)-C2(m)  | 

j^-0  ( (j-m+P2)!i  (j-m+p2-a)b  | 


CQ(0,t)  = 


C* 

1+e(t) 


(1) 


(2) 


*2'  2 

In  Equations  (1)  and  (2),  the  terms  of  the  type  C^k)  are  the  surface 

concentration  of  oxidant  during  the  forward  pulse  of  the  kth  period.  Terms 

of  the  form  C2(k-1)  and  C2(k)  are  the  surface  concentrations  during  the 

reverse  pulses  of  the  (k-1 ) and  kth  periods,  respectively.  Numbering  of 

the  periods  starts  with  j=o  and  C2(-l)  (for  the  first  period)  is  the  bulk 

concentration  C*. 

For  a planar  electrode  and  a simple  reversible  reaction,  the  surface 
concentration  of  oxidant  is  fixed  by  the  potential 


(3) 


4 


where 


e(t)  = exp  j-  [E(t)-E,rJ  j 
and  E,r  is  the  reversible  half  wave  potential. 


Combining  the  potential  wave-forms: 
F ! ( k ) « E,  - kAE-E$w 


during  the  forward  pulse  and 
E2(k)  = E.  - 


k&E+E, 


during  the  reverse  pulse  with  Equation  (3)  giving  the  explicit  surface  con- 
centration to  be  used  in  Equations  (1)  and  (2). 

The  forward  current  can  then  be  written 


Uj) 


oFAdA- 


<t>f  (E$w>  aE»  pr  o) 


and  the  reverse  current 


nFAD^C* 

!r(j)  = ~7^T~  *r  (Esw‘  AE’  p2’  o) 

The  forward  difference  current  is 


Al+(j)  = If(j)  - Ir(j) 


nFAD^C* 


— (*f(j)-*r(j))  - 


nFAD^C 


— Aip+(  j ) 


The  explicit  formulation  of  the  i|»  functions  can  be  written  by  combining 
Equations  (5)  and  (6)  with  Equations  (3)  and  (4)  and  substituting  the  results 
into  Equations  (1)  and  (2). 

Theoretical  square  wave  vol tammograms  were  generated  for  a variety  of 
experimental  parameters.  The  computer  program  explicitly  calculates  and 
stores  the  component  forward  and  reverse  current  functions  and  i^r  rather 
than  merely  the  difference  function  Aip+.  Consideration  of  the  partial  cur- 
rents greatly  facilitates  the  understanding  of  the  results.  The  program 
was  written  in  Fortran  IV  and  was  run  on  Digital  Equipment  Ccrp.  PDP-12 
computer. 


Calculated  square  wave  current-potential  functions  for  several  values 
of  E$w  are  shown  in  Figure  2.  The  difference  current  is  always  an  essentially 
symmetric  peak  shaped  function,  but  for  small  square  wave  amplitude  is  smaller 
than  either  the  forward  or  reverse  currents  (Figure  2A).  The  effect  of  square 
wave  amplitude  on  the  component  currents  and  on  the  difference  current 
clearly  shown  in  Figure  3.  Only  for  n E > ^15  mV  does  the  reverse  current 

bn 

in  the  vicinity  of  the  half-wave  potential  become  negative  and  the  difference 
current  become  larger  in  magnitude  than  the  forward  current. 

The  obvious  conclusion  is  that  analytical  response  can  be  increased  by 
using  larger  square  wave  amplitude,  but,  as  in  pulse  voltammetry,  increased 
amplitudes  lead  to  broadening  of  the  observed  current  peak.  This  effect  is 
shown  in  Figure  4.  The  analytical  choice  must  therefore  be  the  largest  value 
of  square  wave  amplitude  and  step  height  consistent  with  adequate  peak  defini- 
tion and  evaluation.  Of  note  in  Figure  4 is  the  increase  in  the  width  of  the 
current  peak  as  the  square  wave  amplitude  becomes  smaller  than  the  step  height. 

The  step  height  has  little  effect  on  the  height  of  the  current  peak. 

For  a square  wave  amplitude  of  20  mV  the  peak  current  function  Ai|>*ax  goes 
from  0.642  for  an  nAE  = 1 mV  to  0.686  for  an  nAE  = 10  mV  (o  = 5,  p^  = .499, 

P2  s .999).  The  individual  forward  and  reverse  currents,  on  the  other  hand, 
are  markedly  affected  by  the  increase  in  step  height.  This  is  shown  in 
Figure  5.  The  forward  current  function  increases  when  the  step  height  is 
increased  but  the  reverse  current  function  also  increases,  becoming  less 
negative.  The  magnitude  of  change  for  both  currents  is  approximately  the 
same,  so  the  difference  current  remains  relatively  unchanged. 


Figure  5 is  for  constant  x,  therefore  the  sweep  rate  is  increasing  as 
the  step  height  is  increased.  It  is  interesting  to  note  that  increasing  the 
sweep  rate  by  increasing  the  step  height  has  little  effect  on  the  difference 
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current  function.  Increasing  the  sweep  rate  by  decreasing  t will  cause  the 

current  to  increase  as  1//T” for  a reversible  system.  If  we  compare  two  step 

heights  at  constant  sweep  rate,  we  find  a marked  difference  in  Ai|»*  . For 

example,  given  a sweep  rate  of  500  mV/sec  that  we  can  generate  with  a 10  mV 

step  of  20  msec  period  or  a 5 mV  step  of  10  msec  period,  the  normal i?(d  pe.u 
+ -1/2 

current  values  AiL  (t)  ' for  these  are  4.85  and  6.64  respectively, 

max 

The  important  thing  to  remember  is  that  the  peak  current  will  not  only 
depend  on  the  sweep  rate  but  also  on  how  that  sweep  rate  is  generated.  De- 
creasing t to  obtain  a faster  sweep  rate  is  better  than  increasing  AE.  The 
values  of  t that  can  be  used  will,  of  course,  depend  on  the  system  resoonse. 

The  effect  of  measurement  time  is  shown  in  Figure  6.  These  curves  are 
> calculated  for  a symmetric  square  wave  (o=0.5)  and  symmetric  measurement 
times:  p0  = o+p,.  As  the  measurement  time  is  moved  toward  the  beginning  of 
the  pulse,  the  difference  current  function  greatly  increases.  As  Ramaley 
and  Krause  have  shown,  the  increase  in  peak  current  obtained  by  shortening 
the  measurement  time  is  greater  than  the  increase  obtained  by  keeping  the 
measurement  times  at  the  end  of  the  pulse  and  decreasing  the  staircase 
period.  For  the  conditions  of  Figure  6,  the  peak  value  of  the  current 
function  goes  as  p.|"®'®  while  the  peak  current  depends  on  t 

Some  of  the  effects  of  asymmetric  measurement  on  an  asymmetric  square 
wave  are  shown  in  Figure  7.  The  data  for  the  bottom  curve  are  calculated 
for  variable  a with  p1  = o-O.OOl  and  p2  fixed  at  0.999.  For  these  conditions, 
the  peak  current  function  is  a minimum  near  a=0.5  and  increases  essentially 
symmetrically  on  both  sides  of  the  minimum.  The  plot  of  the  maximum  values 
of  a<p+  vs  o becomes  exactly  symmetrical  as  the  ratio  aE/E$w  decreases. 

Figure  8 details  the  individual  forward  and  reverse  current  functions 
for  three  members  of  the  set.  For  o * .1  (Figure  8A)  ^ is  larger  than  Ai|/+ 
because  the  reverse  current  function,  i , has  the  same  sign  as  4/^.  This  is 


h 
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due  to  insufficient  material  being  generated  by  the  forward  pulse  to  allow 
the  reverse  reaction  to  occur  during  the  entire  reverse  pulse.  As  a approaches 
i,  <l'r  gets  more  negative  with  a corresponding  reduction  in  ^ until  at  a = 0.9, 
the  reverse  pulse  is  not  of  sufficient  length  to  regenerate  the  material  con- 
sumed during  the  forward  pulse.  At  this  point  is  at  its  lowest  value  while 
is  at  a maximum  (negatively).  However,  the  change  on  <j>r  is  greater  than 
the  decrease  in  ^ and  thus  the  increase  in  Ai t>+. 

The  data  for  the  top  curve  of  Figure  7 is  again  calculated  for  variable 
o,  this  time  with  fixed  at  0.049  and  = a+0.049.  p-j  and  are  now 
positioned  close  to  the  switching  points  where  < and  iiy  are  maximum  during 
their  individual  steps.  The  peak  function  A^+  is  now  a maximum  at  a = .5, 
decreasing  almost  symmetrically  on  each  side.  Again  the  plot  becomes  exactly 
symmetrical  at  large  square  wave  amplitude  to  step  height  ratios. 

Individual  forward  and  reverse  currents  for  three  members  of  this  set 
are  shown  in  Figure  9.  The  trends  for  the  individual  current  function 
and  ^ are  the  same  here  as  in  Figure  8.  ^ decreases  and  i)»r  increases 

(negatively)  as  o approaches  0.9.  However,  in  contrast  to  Figure  8 where 
the  negative  magnitude  of  i|*r  changes  dramatically  as  o goes  from  0.1  to 
0.9,  the  change  in  i|<r  here  is  relatively  small. 

Figure  It) illustrates  the  information  that  can  be  gained  by  a study  of 
the  individual  forward  and  reverse  current  functions.  This  is  a vertical 
comparison  of  two  members  in  Figure  7.  Here  a and  p.j  are  the  same  for  both 
figures,  equal  to  .1  and  .099  respectively  and  p ^ is  varied  from  close  to 
the  switching  point  (0.149)  Fig.l0Ato  near  the  end  of  the  step  (0.999)  Fig.  10B. 
From  examination  of  the  reverse  current  function,  it  is  obvious  that  the 
current  during  the  reverse  step  has  actually  changed  direction.  Directly 
after  the  forward  step  4>r  is  opposite  in  sign  to  while  near  the  end  of 
the  step  i i>f  and  has  the  same  sign.  This  indicates  that  the  material 
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generated  by  the  forward  pulse  has  been  consumed  during  the  first  part  of 

the  reverse  puise  and  in  order  for  the  Nernst  equation  and  Fick's  laws  of 

diffusion  to  be  obeyed,  the  current  near  the  end  of  the  step  must  now  flow 
in  the  same  direction  as  4^. 

There  is  no  obvious  reason  to  output  the  difference  current  in  which 

the  current  at  during  the  jth  cycle  is  differenced  with  the  current  at 

P2  during  the  same  staircase  period.  Since  we  have  calculated  and  stored 
the  partial  current  components,  it  is  a simple  matter  to  examine  other  dif- 
ferencing schemes,  such  as  the  backward  difference,  Ai|T(j)  = j ) - 4y(j-l). 
This  difference  function  has  somewhat  larger  maximum  values,  especially  at 
smaller  square  wave  amplitudes,  as  can  be  inferred  from  an  inspection  of  the 
potential  dependence  of  the  partial  components  shown  in  Figure  2.  The  apparent 
peak  potential  is  shifted  by  approximately  one-half  the  step  height,  relative 
to  the  maximum  of  the  forward  difference.  Some  representative  and  comparative 
values  of  the  maximum  peak  current  function  U41  1,  peak  potentials  and  peak 

f fTlaX 

widths  for  forward  and  backward  differences  are  shown  in  Table  I. 

The  similarity  of  this  technique  to  pulse  polarography  is  obvious, 
although  as  has  been  demonstrated,  variations  in  points  of  measurement 
affect  the  shape  of  the  wave.  Nevertheless,  the  equations  developed  result 
in  wave  shapes  identical  to  those  for  differential  pulse  polarography  (1) 
at  slow  scan  rates  and  measurements  parameters  identical  to  those  employed 
in  pulse  polarography. 

Work  Is  in  progress  to  develop  the  theory  to  include  the  case  of  slow 
electron  transfer  processes,  and  first-order  chemical  processes,  including 
preceeding,  following  and  catalytic  reactions. 

Experimental  verification  of  the  theory  presented  here  and  demonstration 
of  the  analytical  utility  of  the  method  at  the  dropping  mercury  electrode 
will  be  forthcoming  (6). 

I - — — - 
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Figure  1 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


FIGURE  CAPTIONS 


Square  wave  waveform  with  associated  nomenclature,  t is  the 
step  width  for  the  base  staircase. 

Individual  and  difference  current  functions  for  three  square 
wave  amplitudes.  nAE  = 5 mV,  o = 0.5,  p1  = 0.499,  r2  = 0.999 

(A)  nE$w  = 5 mV 

(B)  nE$w  = 15  mV 

(C)  nE$w  = 30  mV 

Solid  line  is  difference  current  function,  Atj,+ 

Individual  and  difference  current  functions  at  the  peak 
(a^x)  vs  square  wave  amplitude.  nAE  = 5 mV,  o = 0.5, 
p-|  = 0.499,  p2  = 0.999. 

( ■ ) Aip 

( • ) at  peak 
(A)  ij)r  at  peak 

Peak  width  at  half  height  vs  square  wave  amplitude.  nAE  = 

5 mV,  a = 0.5,  p-j  = 0.499,  p2  = 0.999. 

Individual  forward  and  reverse  current  functions  as  a function 
of  AE.  nEgw  = 20  mV,  a = 0.5,  p.|  = .499,  p2  = .999 

(A)  nAE  = 10  mV 

(B)  nAE  = 5 mV 

(C)  nAE  = 2 mV 

(D)  nAE  = 1 mV 

Effect  of  measurement  time  on  A*+.  nAE  * 5 mV,  nE$w  = 20  mV, 
a = 0.5,  p2  * 0 + pi 


(B)  P]  = 0.1 

(C)  Pl  = 0.3 

(D)  P]  = 0.499 


f 
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Figure  7. 


Figure  8. 


Figure  9. 


Figure  10. 


Peak  current  function,  a/.  vs  o.  nAE  = 5 mV,  nE  =20  mV. 

max  SW 

(♦)  P]  = a - 0.001,  P2  = 0.999 

(A)  P]  = 0.049,  p2  = o + 0.049 

Effect  of  asymmetric  measurements  on  asymmetric  square  waves. 
nAE  = 5 mV,  nEsw  = 20  mV,  P.j  = o-.OOl,  p2  = .')99. 

(A)  o = 0.1 

(B)  o = 0.5 

(C)  o = 0.9 

Symmetry  of  waveform  shown  in  inset,  (x)  indicates  current 
measurement. 

Effect  of  asymmetric  measurements  on  asymmetric  square  waves. 
p.|  = 0.049,  p2  = o + 0.049,  other  conditions  and  legend  same 
as  Figure  8. 

Change  in  current  function  with  sampling  time  nAE  = 5 mV, 
nE  = 20  mV,  a =0.1,  p,  = 0.099 

bW  ! 

(A)  p2  = 0.149 

(B)  p2  = 0.999 
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